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Abstract 

We present a model for magnetic structured coronae 
above accretion discs. On the shortest timescales, spatially 
and temporally correlated coronal flares can explain X- 
ray temporal and spectral variability observed in Seyfert 
galaxies. In particular, power density spectra, flux-spectral 
index and flux- variance correlations are naturally accoun- 
ted for by the model. More dramatic spectral variations 
(i.e. state transitions in GBHC) are associated with pa- 
rameters varying on longer timescales, such as accretion 
rate, coronal strength or geometry of the inner disc. In 
the framework of the standard Shakura-Sunyaev accretion 
disc theory, here we discuss why energetically dominant 
coronae at low accretion rates are ideal sites for launch- 
ing powerful MHD driven outflows. Then, if the outflow is 
radiatively inefficient, then so is the source overall, even 
without advection being relevant for the dynamics of the 
accretion flow. This could be an alternative scenario for 
LLAGN and GBHC in their low/hard state, and may have 
consequences for our understanding of the accretion his- 
tory of the universe. 
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1. Introduction 

The hard X-ray spectra, the properties of the reflection 
features and of rapid time variability in accreting black 
holes can all be considered as indications of the presence 
of a structured, hot, optically thin component in the inner 
part of the accretion flow (the corona), situated above a 
cold, geometrically thin, disc. Many models of the accre- 
tion disc corona have been proposed in recent years, which 
are able to fit the observed time averaged spectra in terms 
of Comptonization of soft photons in the hot corona. How- 
ever, many uncertainties regarding the actual geometry of 
the inner accretion flow remain unsolved (see Done 2001 
for a review). Here we present a model for a structured, 
magnetic corona above a standard, geometrically thin and 
optically thick disc, that can satisfy observational con- 
straints on both long and short time variability of these 
systems. 



Let us briefly sketch the overall energetics of an accre- 
tion disc-corona system. The total power released by the 
accreting gas is defined as L = iriL^dd — 47rGMm p mc/<7T. 
To be fairly general, we assume that the coronal power 
generated by the disc (which is a fraction / of the total: 
L c = /fh-ZjEdd) can be either dissipated locally to heat 
the corona, and ultimately radiated away as hard X-rays 
with a luminosity Ljj = (1 — v)Lc, or used to launch an 
outflow with power Lj = r)L c . The main characteristic of 
our model are the following: in the hot corona the energy 
is stored in a strong, highly intermittent magnetic field 
(amplified in the turbulent disc and buoyantly expelled 
in the vertical direction) . Magnetic energy dissipation oc- 
curs at the smallest end of the turbulent-energy cascade. 
Such small flares heat the corona (with a power Lh) and 
can trigger an avalanche in their immediate neighborhood 
(Poutanen & Fabian 1999), creating a bigger active re- 



gion, and producing the observed flares in the lightcurves. 
In Sec. [2[ a stochastic model (the so-called thundercloud 
model, Merloni & Fabian, 2001, MF1) of the short-time 
variability is discussed. 

On longer timescales, variations of the fraction of power 
released into the corona, /, and of the fraction of coronal 
power used to launch an outflow, rj, may be associated, 
for example, to changes of the accretion rate and/or of 
the geometry of the inner disc. In Sec. ^ we show that, 
under reasonable assumptions on the nature of the disc 
viscosity, the strength of the corona / increases as the 
accretion rate decreases. Also we discuss reasons why an 
energetically dominant corona can be the site where pow- 
erful outflows are produced. 



2. Short timescales: magnetic flares 

Seyfert 1 galaxies show strong temporal variability in the 
X-ray band, with a Power Density Spectrum (PDS) char- 
acterized by a "red noise" power-law: P(f) — Af~ a with 
a ~ 1.5 (see e.g. Nandra, 2001). The dynamical timescale 
of such systems, t dyn = 9 x lO 3 r 3 / 2 (M B H/lO 8 M ) s, is 
such that good time-resolved spectroscopy can directly 
provide information about the heating mechanism of the 
corona. In fact, from the observational point of view there 
is a growing evidence that during a flare the X-ray spec- 
trum becomes softer as the (2-10 keV) flux increases (see 
e.g. Vaughan & Edelson 2001). Also, a remarkable linear 
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correlation between variability and flux on all timescales 
has recently been reported (Uttley & Mc Hardy 2001). 

The coronal magnetic field intensity depends on the 
dissipation rate and on the field geometry. Here we as- 
sume, on the basis of the strong spectral and temporal 
variability, that the corona is patchy and the magnetic 
field concentrated in a number of small active regions 
where dissipation occurs ( Haardt, Maraschi fc Ghiscllin] 
1994|; |Sjtcrn et al. 1995|; pi Matteo, Celotti fc Fabian 1999[)' 



We assume that the luminosity of an active region is 
determined by the size of the avalanche, and therefore 
scales as the size of the region: L(R a ) = C\R® . Due to the 
fractal nature of reconnection and avalanche generation, 
we also assume that the number of active regions of size 
R a at any time, is given by dN(R a )/dR a — C^Ra 1 '■ The 
PDS slope, a, is then related to the indices D and p by 
the relation p = 2D + 3 — a. The constants C\ and C2 can 
be fixed imposing the overall normalization for the total 
corona covering fraction and the average hard luminosity 
(see MF1 for details). The covering fraction C is related 
to the total variability level as illustrated in the example 
in Fig. | 
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Figure 1. Average PDS of ten segments of a simulated 
lightcurve (each 3 x 10 seconds long), for fixed coronal 
optical depth, r = 1.5 and D — 1. From bottom to top 
the different curves correspond to a covering fraction C = 
0.03 (orange, open circles), 0.01 (purple, open squares), 
0.005 (cyan, open triangles), 0.002 (blue, filled stars), 
0.001 (green, filled circles), 0.0003 (red, filled squares) and 
0.0001 (black, filled triangles) . On the right, close to each 
curve, the corresponding fractional rms is reported. 



As discussed in MF1, we have simulated AGN X-ray 
lightcurves fixing a = 1.5 (in order to reproduced the 
observed properties of the PDS), and studied the time 
and spectral variability properties of the model. 

We assume that the X-ray spectrum is produced by 
thermal Comptonization in each spherical active region 
of size R a , lifted above the disc at a height H a and with 
Thomson optical depth r. The active regions illuminate 



the cold disc: part of the flux is reflected and part is ab- 
sorbed and reprocessed. We assume that most of the ac- 
cretion power is dissipated in the corona (/ > 0.7), so 
that reprocessed radiation is the main source of soft pho- 
tons in the active regions. For large and luminous active 
regions the amount of soft radiation intercepted by an ac- 
tive region increases, and the spectrum is softer. This is 
illustrated in the example of Fig. |2|. 

We calculate the spectrum using an analytic approxi- 
mation for thermal Comptonization (see MF2), fixing H a = 
i?s = 2GMbh/c2, and we study the correlations between 
spectral properties and total luminosity (in the 3-10 keV 
band). We have fixed Mbh = 10 7 Mq and an accretion 
rate of 10% of the Eddington one. In MF1 we show in 
detail how such a model can reproduce the observed lu- 
minosity - spectral index correlation. In the following sec- 
tion, instead, we will focus on the luminosity-variability 
correlation (see also Merloni 2001). 

2.1. The luminosity- variability correlation 



Recently, Uttley fc McHardy (2001) have discovered a 
remarkable linear correlation, on all timescales (Uttley, 
private communication), between variability and flux in 
many compact accreting systems. 

We have analysed one of our simulated lightcurves in 
order to look for a similar correlation. With the parame- 
ters D and C fixed to the values of 1 and 0.001, respec- 
tively, we have produced a lightcurve of the total duration 
of about 3 x 10 6 s, with resolution of 1 minute. 

Then, keeping fixed the time resolution, we calculated 
mean and variance for segments of different length. The 
results, binned in luminosity intervals, are shown in Fig. ^, 
where the length of the segments T is measured in min- 
utes. The variance corresponds to the integral of the PDS 
shown in Fig. [l] from frequencies 1/(T x 60) to 1/60 (Hz). 

A clear correlation is evident in all cases, consistent 
with a linear one with an offset. |Uttley fc McHardy (2001 \ 
argues that the correlation rules out any possible shot 
noise model. Still, we are able to reproduce it within our 
thundercloud model, which is indeed a generalized shot 
noise model. One possible explanation for this lies in the 
way the lightcurve is 'randomized'. Once the average cov- 
ering fraction C and the distributions indexes D and p 
are fixed, the total expected number of active region at 
any time N tot (t) is calculated, then the actual (random- 
ized) number, A/"tot(i)> is picked from a poissonian distri- 
bution that has N tot (t) as its mean. Then, starting from 
the smallest flares, we calculate in similar fashion the num- 
ber of active regions of every size in such a way as to ensure 
that their total is exactly A/tot (£)■ Therefore, random vari- 
ations in A/tot (t) (which correspond to random variations 
in the covering fraction C) are transferred to flares of all 
sizes, and of all timescales, thus enforcing a correlation be- 
tween the overall flux level and the amount of variability 
on all time-scales. This merely reflects the fact that large 
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Figure 2. The inner corona is mapped into a square. For simplicity, we chose D — 2, so that the pixels have equal 
luminosity, and the total luminosity scales with the covered area. The region of the corona active at any time are 
represented by filled squares, whose size is distributed as a power-law, with index p — 2. Each region is assumed to 
have the same height above the disc (H a = Rg) and optical depth (r = 1). The color code shows the spectral index of 
the filled active regions, calculated self- consistently taking into account thermal Comptonization of disc photons (both 
intrinsic and reprocessed) and of synchrotron emission in the active regions themselves (see MF1). It ranges from 
r = 1.62 (pink) to T = 2.21 (red). Clearly the more luminous state (right hand side, for which the instantaneous 
covering fraction is larger), is dominated by active regions with soft spectra. 
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Figure 3. Luminosity dependence of the variance calculated 
for segments of a simulated lightcurve (D — 1, C — 10 -3 
and t = 1.5) with a fixed temporal resolution of 60 sec- 
onds. The duration of the different segments ( in minutes ) 
is written in each panel. The fitting functions are of the 
form a — fc(£3_io — La), with k and Lq constants. 



active regions are indeed made of a collection of correlated 
microflares (avalanches) . 

3. Long timescales: low-luminosity systems 

On timescales much longer than the dynamical one, the 
system evolution is regulated by the global properties of 



the accretion flow. Here we want to discuss the implication 
of our coronal model for the long-term properties of sys- 
tems accreting at low rates (as compared to the Eddington 
one). In particular we will consider Galactic Black Hole 
Candidates (GBHC) in their low/hard state and the so- 
called Low Luminosity AGN (LLAGN). These two classes 
of systems have many common features: GBHC have very 
weak (or absent) quasi-thermal spectral component in the 
soft X-ray/EUV spectral range; hard (photon index L = 
1.5 — 2) X-ray power-law spectra, rolling over at ~ 100 keV 
(evidence of thermal Comptonization in a optically thin 
medium); compact (unresolved) flat or inverted radio core, 
probably extending to the NIR or optical regime (??) and 
a clear temporal correlation between the radio and hard 
X-ray fluxes. LLAGN, in turn, show a remarkable weak- 
ness of the Big Blue Bump ( Ho 1999 ); hard X-ray spec- 
tra (r = 1.6 — 1.8); compact radio cores nearly ubiqui- 
tously, often accompanied by 'jet-like' features; radio spec- 
tra with turnover frequencies of the order ~ 0.1 — 100 GHz 
and radio loudness anti-correlated with accretion rate (for 
nearby ellipticals with reliable mass estimates; Ho 2002). 



All these observational evidences suggest that outflows 
are common features of low m systems with hard X-ray 
spectra. 
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3.1. Corona dominated accretion discs 

As shown in Merloni & Fabian (2002, MF2), if angular mo- 
mentum transport in a standard, geometrically thin and 
optically thick accretion disc is due to magnetic turbu- 
lent stresses, the magnetic energy density and the effec- 
tive viscous stresses inside the disc are proportional to the 
geometric mean of the total (gas plus radiation) and gas 
pressure: 

— = Pmag — CtQ J -Ptot-Pgas, (1) 
07T 

with ao constant (of the order of unity) . 

Then, the fraction of gravitational power / released 
in a magnetic corona increases as the accretion rate de- 
creases, because the disc is more and more gas pressure 
dominated even in its inner parts. When the disc is com- 
pletely gas pressure dominated (at accretion rates smaller 
than the critical value rh cr i t ~ 0.1(o:oMbh/-M©) _1 / 8 )) the 
strength of the corona depends on the poorly understood 
mechanisms of vertical flux tube transport in the disc, but 
it is likely that / > 0.7 (see MF2 for a more thoroughly 
discussion) . 



3.2. MHD JETS FROM MAGNETIC CORONAE 

The magnetic energy density in the corona is given by 
IP _ _3i 
8tt " 



4irRlN tot c \v diss J ' (2) 
where «diss is the dissipation velocity and depends on the 
uncertain nature of the reconnection process (MF1), and 
can be assumed to be of the order Vdiss ~ 0.01c. 



Models and s imulations o f jet production ( Blandford 
fc Znajck 1977| ; |Meier 1999| ) show that it is the poloidal 
component of the magnetic field which mainly drives the 
production of powerful jets, and the output power can be 
expressed as (Livio, Ogilvie & Pringle 1999) 
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Here R cor is the size of the region where most of the coro- 
nal power is dissipated and Q is the typical angular veloc- 
ity of the magnetic field lines. 

If H a is the typical coronal flux tube scaleheight (height 
of a reconnection site), we have B p /B ~ H a /R cor , and 
therefore 
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v ^diss / \ -*£cor / \ C 

which in turn gives, for the fraction of coronal power that 
goes into the MHD outflow 

/ . . 2 / . \ -i 

2 / ^diss \ / Rcnr \ ( C 

V ' 
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(5) 



By inspection of Eq. (g) we can conclude that the outflow 
power is stronger if the coronal scaleheight is large with 
respect to the distance from the central source. This would 



help in increasing the relative strength of the poloidal com- 
ponent of the magnetic field, that is the one ultimately 
responsible for the powering of the jet. 

As an example, for Vdiss — 0.01c, R COI ~ 7Rs, H a ~ 
2Rs and ft = Qk(Rcoi), we obtain ?7mhd — 0.55: the MHD 
jet can carry away a substantial fraction of the coronal 
power. 

4. Conclusions 

We have presented a model to explain spectral and tem- 
poral variability on the smallest timescales in the X-ray 
emission from Seyfert Galaxies and GBHC. We have simu- 
lated X-ray light-curves that reproduce the observed PDS 
properties and the spectral variability. The basic geomet- 
ric properties of the corona we propose are the following: 

— The corona must not be uniform, but structured and 
heated intermittently (flares); 

— The spatial and temporal distribution of the flares are 
not random, but proceed in correlated trains of events 
(avalanches) ; 

— The size of the avalanches determines the size of the 
active regions, which are distributed as a power-law; 
larger avalanches are more luminous and have softer 
spectra. 

On longer timescales, the evolution of the corona is 
governed by the evolution of the accretion rate and/or 
of the inner disc geometry. At low accretion rate, the 
strength of a magnetic corona produced by buoyant mag- 
netic flux tube amplified in an underlying standard ac- 
cretion disc increases. If the energy in the corona, as we 
suggest, is stored in the magnetic field, and the height of 
a reconnection site is much larger than its size, which is 
of the order of the disc thickness, powerful MHD outflows 
can be launched from the inner corona. 
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